SCIENTIFIC 

REPORTS 



s 




OPEN 



SUBJECT AREAS: 
GEOCHEMISTRY 
PETROLOGY 

Received 
21 January 2014 

Accepted 
14 May 2014 

Published 
17 June 2014 



Correspondence and 
requests for materials 
should be addressed to 
S.C.H. (huangl7@fas. 
harvard.edu; huangs@ 
alum.mit.edu) 



Missing Lead and High ^He/'^He in 
Ancient Sulfides Associated with 
Continental Crust Formation 



Shichun Huang\ Cin-Ty A. Lee^ &Qing-Zhu Yin^ 



^ Department of Earth and Planetary Sciences, Harvard University, ■^Department of Earth Science, Rice University, ^Department of 
Earth and Planetary Sciences, University of California, Davis. 

Major terrestrial reservoirs have Pb isotopes more radiogenic than the bulk silicate Earth. This requires a 
missing unradiogenic Pb reservoir, which has been argued to reside in the lower continental crust or 
dissolved in the core. Chalcophile element studies indicate that continent formation requires the formation 
of sulfide-bearing mafic cumulates in arcs. Because Pb, but not U, partitions into sulfides, we show that 
continent formation must have simultaneously generated time-integrated unradiogenic Pb reservoirs 
composed of sulfide-bearing cumulates, now recycled back into the mantle or stored deep in the continental 
lithosphere. The generation of such cumulates could also lead to coupled He-Pb isotopic systematics because 
^He is also produced during U-Th-Pb decay. Here, we show that He maybe soluble in sulfide melts, such that 
sulfide-bearing cumulates would be enriched in both Pb and He relative to U and Th, "freezing" in He and 
Pb isotopes of the ambient mantle at the time of sulfide formation. This implies that ancient sulfide-bearing 
cumulates would be characterized by unradiogenic Pb and He isotopes (high-^He/^He). These primitive 
signatures are usually attributed to primordial, undifferentiated mantle, but in this case, they are the very 
imprint of mantle differentiation via continent formation. 



The geochemical behavior of Pb in the Earth is not fully understood. Most upper crustal and upper mantle 
reservoirs have Pb isotopes more radiogenic than bulk Earth estimates, manifested by silicate reservoirs 
plotting to the right of the geochron on a 206p|^/204p|^ versus ^°^Pb/^°^Pb diagram. This imbalance, known as 
the "Pb paradox", means that a missing reservoir with unradiogenic Pb is required to balance the radiogenic Pb\ 
One solution is that the missing unradiogenic Pb resides in the present-day lower continental crust^. However, 
most lower crustal xenoliths are not sufficiently unradiogenic to satisfy this mass balance^'^'^. Another hypothesis, 
motivated by experimental studies, is that Pb dissolved into the Earth's core early in Earth's history^'^. If Pb indeed 
segregated into the core, the U-Pb systematics of the bulk silicate Earth would indicate that the core formed 80- 
140 Myr after planetary accretion^'^, but this core formation age is much later than the < 30 Myr after accretion 
based on ^^^Hf-^^^W isotopic systematics, which is uniquely sensitive to core formation^'^ (unlike the U-Pb system, 
which can be fractionated by processes unrelated to core formation). Alternatively, this young U-Pb age of the 
bulk silicate Earth, 80-140 Myr after Earth's accretion^'^, has been used to argue for a late accretion of volatile 
elements, such as Pb, to the Earth^°'^\ Due to these inconsistencies, additional hypotheses should be explored. 

Sulfide minerals are another potentially important carrier of Pb^^. A recent study by Lee et al.^^ showed that 
pyroxenite cumulates formed during the differentiation of arc magmas are enriched in sulfides and therefore 
enriched in strongly chalcophile (sulfide-loving) elements like Cu. The bulk continental crust is felsic in com- 
position because a large amount of mafic cumulates were formed and then removed from the continental crust by 
delamination^^'^^. Lee et al.^^ showed that the Cu content of the continental crust is depleted relative to basaltic 
magmas and even the mantle, which implies that the complementary mafic cumulates should be sulfide-bearing. 

The amount of S, in the form of sulfide, which must be missing from the continental crust, can be evaluated by 
considering the compositional effects on sulfide solubility during magmatic differentiation. Assuming the initial 
building blocks of continents are basalts and that almost all basalts are sulfide- saturated in their mantle source 
regions, the initial S contents of these juvenile basalts will range from 1600-2000 ppm^^. Because sulfide solubility 
decreases substantially with increasing Si02 and decreasing FeO contents, by the time the basalts evolve to 
andesitic compositions typical of bulk continental crust, the S content at sulfide solubility has fallen to 
<200 ppm^^ (solubility may rise with further magmatic evolution if oxygen fugacities rise, but it is the andesitic 
compositions that are of interest here as such compositions most closely match the average composition of the 
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continental crust). Based on the foregoing analysis, possibly —90% of 
the original S must be either sequestered in the form of sulfide- 
bearing cumulates or lost during volcanic degassing^^'^^, leaving 
behind a S-depleted continental crust^^. 

The presence of a sulfide-bearing mafic reservoir, a required by- 
product of continent formation, has important implications for Pb. 
Pb is a moderately chalcophile element. Depending on the Pb par- 
tition coefficient in sulfide (14-40)^^, these sulfide-rich cumulates 
may sequester 5-15% of the parental magma's Pb, the rest remaining 
in the continental crust^^. Because these sulfide-bearing cumulates 
sequester negligible amounts of U and Th, they would be expected to 
"freeze" in the Pb isotopic composition of the mantle at the time of 
their formation, generating a reservoir for unradiogenic Pb. In fact, 
unradiogenic Pb isotopic compositions, plotting on the left side of 
the geochron in a ^°'Pb/'°^Pb vs. '°'Pb/''^Pb diagram, have been 
found in sulfide inclusions from abyssal peridotites in the Atlantic 
Ocean^^. 

To explore this hypothesis further, we turn to He because ^He is 
also a decay product of U-Th-Pb alpha-decay systems. Sulfides 
formed during magmatic differentiation are likely to be initially in 
the liquid state due to their low melting temperatures^^. The question 
arises as to whether noble gases like He might be soluble in sulfide 
liquids. As long as there is a finite solubility of He in sulfide liquids, 
regardless of how small, and because U and Th do not partition at all 
into sulfides, we hypothesize that sulfides could effectively "freeze" in 
the He isotopic compositions of the mantle at the time of formation. 
If formed early in the Earth's history, these sulfides would be char- 
acterized by low 2°^Pb/'°^Pb and high ^^^Ph^nh^ (the ratio of 
radiogenic ingrowth of ^°^Pb and ^^^Pb since the formation of the 
Earth2\ and defined as ^^'PhT'Vh'' = (^^^Pb/^^^Pb- 10.297)/ 
C°'Pb/'°^Pb-9.307)) and 'He/'He, thus Pb-He isotope systematics 
should be coupled if sulfide control is important, generating negative 
206p|3/204p|^_3Hg/4p^g and positive ^^^Pb^/'^'Pb^-'He/^He correlations. 
It is thus worth evaluating whether He is soluble in sulfide liquids. 

Noble gas studies of sulfides provide some insight. Using step- 
heated degassing experiments, Yang and Anders^^ showed that 30- 
70% of the total noble gas budget of FeS sulfides is housed in the high 
temperature fraction, suggesting that noble gases can dissolve directly 
into the sulfide structure. Given that sulfide liquids are less dense and 
hence likely to contain more atom- scale void space than sulfide 
minerals, the fact that small amounts of noble gases can dissolve into 
the structure of sulfide minerals suggests that noble gases could be 
even more soluble in liquid sulfides than sulfide minerals. 

The solubility of noble gases in basaltic melts has been experiment- 
ally investigated^^'^^, but as far as we know, there are no solubility 
experiments of noble gases in sulfide liquids. For this reason, we 
attempt here to infer He solubilities indirectly. Carroll and 
Stolper^^ showed experimentally that the noble gas solubilities in 
silicate melts correlate well with ionic porosity, allowing them to 
predict solubilities for those noble gases that have not been studied 
in detail. Unfortunately, we are unable to apply this method to sulfide 
melts because a completely different and unknown parameterization 
is needed given the stronger covalent bonding seen in sulfides than in 
the metal sites of silicate melts. We thus adopt a different approach 
based on the work of Uhlig^^, Blander et al.^^ and Lux^^. Assuming 
that He, being a non-polarizable noble gas, does not interact with the 
liquid solvent, the free energy change associated with dissolution 
should be related to the surface energy of the void created by the 
gas^^'^^. Taking the macroscopic surface tension y (mj/m^) as an 
approximation of the surface energy associated with the void, the 
solubility of a noble gas in a liquid has been shown semi-empirically 



approximation of the size of the void. Consequently, the partitioning 
between sulfide and basaltic melts has the following relationship: 
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where Csf and Csu are 

the noble gas concentrations in sulfide and basaltic melts and y^y 
and Ysii are the macroscopic surface tensions of sulfide and basaltic 
melts. Using the radius of a He atom (1.08 A) in six-fold coordina- 
tion^^ and surface tensions of —370 mj/m^ and —550 mj/m^ for 

basaltic and sulfide melts at T = 1300°C'°''^ D^^'^' is estimated to 
be —0.3. It has been shown that if small amounts of O are dissolved in 
the sulfide melt, as would be expected at the oxygen fugacity condi- 
tions typical of the uppermost mantle (near the quartz-fayalite-mag- 
netite buffer), the wetting angle and by implication the surface 
tension of a sulfide melt decreases considerably^^. This would 
decrease the difference between sulfide and basaltic melt surface 

tensions, driving D^^i^^^ towards greater values. 

Because He is almost perfectly incompatible in silicate minerals^ \ 
most of the He budget within a partially molten system will be con- 
tained in the melt phases. Thus, a D^^^^ value of 0.3 or higher impUes 
that non-trivial amounts of He can dissolve in sulfide melts. If sulfide 
melts, owing to their high densities compared to basaltic melts, can 
be efficiently segregated from basaltic melts, these sulfide melts 
become reservoirs characterized by low (U + Th)/He and (U + 
Th)/Pb because neither U nor Th partitions into sulfides^^. Such 
sulfide melts can be sealed as inclusions in pyroxenes and olivines 
that form the mafic cumulates. Since He is insoluble in the host 
silicate mineraP^ once the sulfide liquid is trapped. He is unlikely 
to escape the inclusion, even if sulfide liquids crystallize into sulfide 
minerals and temperatures remain high enough for He to diffuse 
rapidly (akin to silicate minerals being entombed as inclusions in 
diamond and preserving their isotopic signatures). In summary, we 
speculate that sulfide inclusions in mafic cumulates freeze in the 
'Heme, ^«^Pb/'°^Pb and 207pb*/206pb* ^^^^le ambient mantle at the 
time of formation. This has important implications for the isotopic 
evolution of these cumulates. 

We now model the secular evolution of '°'Pb/'°^Pb, ^^'^hT'Fh'' 
and ^He/^He in the mantle for closed system and continuous degass- 
ing (Fig. 1). As an upper bound on the effect of He degassing, we only 
consider constant He degassing rate scenarios. During degassing, only 
He, but not U, Th or Pb, is allowed to leave the system (e.g., into the 
atmosphere) because U, Th and Pb transported to the crust are then 
recycled to the mantle, and thus remain within the condensed part of 
the Earth. The initial 'He/^He is taken as 120 R/Ra'^ so that the 
hypothetical present-day primitive mantle has a ^He/^He of 54 R/ 
Ra (Fig. Ic). The co-evolution of He and Pb isotopic compositions 
should thus form a negative ^^^Pb/^^^^Pb-^He^He trend and a positive 
'°'Pb*/'°'Pb*-'He/^He trend, the steepness of the slope depending on 
the integrated extent of degassing (Figs. 1, 2). Losing 98.5% of the 
primordial ^He from the Earth's interior results in the present-day 
MORB 'He/^He value of 8 R/Ra. For comparison, the highest 'He/'He 
ever found on the Earth (in Baffin Island picrites) is 50 R/Ra^^ which 
requires almost no degassing or that the He isotopic composition of 
the mantle was frozen in billions of years ago (Fig. 2). Mathematically, 
our model for ^He/^He evolution of partially degassed primitive man- 
tle is similar to those discussed elsewhere^^"^^, but our model is the 
first that may provide a clear link between He-Pb isotopic evolution 
and the differentiation of continental crust. 

We compiled Pb and He isotopic data from 16 oceanc island basalt 
(OIB) and continental flood basalt (CFB) localities in Fig. 2 and Table 
S2. Within individual OIB-CFB suites. He and Pb isotopes appear to 
be decoupled, most likely because OIB-CFB magmas are mixtures 
between endmembers having vastly different He/Pb ratios due to 
previous or contemporaneous degassing events^^. The global OIB- 
CFB data display a boot-shape distribution in a 206p|^/204p|^ 
^He/^He plot, but the average of each OIB-CFB suite shows negative 
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Figure 1 | Time evolution of the primitive mantle and a series of partially 
degassed primitive mantle reservoirs. For the primitive mantle, the 
following initial values (at 4.567 Ga, all atomic ratios) are used: ^^^U/^°^Pb 
= 16.96, ^^^Th/^^^Pb = 41.29, ^/^"^Pb = 0.0215, %me = 1.66 X 10-^ 
and 206p|^/204p|^ ^ 9.307, so that at present-day, the primitive mantle has 
238u/204p|^ = 8.35 and'He/^He = 54 R/Ra. For partially degassed primitive 
mantle reservoirs, only He, but not U, Th or Pb, is allowed to leave the 
system. We use a continuous degassing process, in which the rate of ^He 
loss is constant at each case and ^He/^He is not fractionated during 
degassing. 207p|^>i</206p|^* ^j^^ ratio of radiogenic ingrowth of ^°^Pb and 
2°^Pb since the formation of the Earth'^ and is defined as 207pb*/206p|^* = 
(207p|^/204p|3_^Q 297)/(206p^/204p|^_9 3Q7) r^^^ pcrccutages labeled in Panel 
c are the present-day ^He budget left in each partially degassed primitive 
mantle reservoir, and the 100% model line represents the undegassed 
primitive mantle. The calculation was done step -wisely with a step of 
10 Ma, and the details are in Table S 1 . The two horizontal trends in Panel c 
label the present-day average MORB ^He/^He value (8 R/Ra) and the 
highest 'He/^He observed (50 R/Ra). 



'°'Pb/'°^Pb-'He/^He and positive 2°^Pb*/'°'Pb*-'He/^He trends, with 
the high ^He/^He associated with ancient Pb isotopic signature, i.e., 
low ^^'Fhr^Vh and high ^^'^h'^ r^^h"" (Figs. 1, 2b, 2c). We note that 
due to crustal contamination or post-eruptive radiogenic ingrowth of 
^He, the average ^He/^He for each OIB-CFB suite likely represents a 
minimum bound on the isotopic signature of the primary magma. 
Thus, a more informative approach is to plot the maximum observed 
^He/^He at each OIB-CFB suite. As shown in Figs. 2d and 2e, max- 
imum observed values also form negative ^^^^h/^^^I'h- ^He/^He and 
positive 207p|^>k/206p|^H«_ ^Hq/^iiq trends, implying that such correla- 
tions are mantle signatures. 



High ^He/^He is usually explained as a signature of undegassed and 
unprocessed mantle. For example, Jackson et al.^^ observed that high 
^He/^He is associated with ancient Pb isotopic signature in Baffin 
Island and West Greenland lavas, and they proposed that this reflects 
sampling of undegassed primitive mantie. Recently, alternative views 
have entertained the possibiUty that high ^He/^He might instead 
reflect ancient residues of mantle melting owing to preferential parti- 
tioning of U relative to He into the melt^\ but it is not clear how such 
a scenario would affect time-integrated Pb isotopes. Sulflde-rich 
cumulates formed during the formation of continental crust early 
in the Earth's history^^ would be characterized by low 206p|^/204p|^ 
and high 'He/'He and ^^'vhT'^h'^ as discussed above (Figs. 1, 2). 
The great majority of these cumulates are thought to have delami- 
nated back into the mantie or been stored as small pockets within the 
continental lithosphere^l If these cumulates are tapped by plumes, by 
small scale convective instabilities or are preferentially melted during 
lithospheric heating/extension, the coupled Pb-He isotopic system- 
atics would emerge in OIB systematics as can be seen in Figure 2. 
These signatures would not be expected to appear in mid-ocean ridge 
basalts because the higher extents of melting beneath ridges would 
dilute the pyroxenite signature by the signature of ambient mantie. 

If ancient sulflde-bearing cumulates exist in the convecting mantle 
or continental lithospheric mantle, we might expect to see coupled 
Pb-He signatures in continental intraplate magmas as well. One of 
the few continental intraplate localities that have been extensively 
studied in terms of Pb and He isotopic systematics is the 
Yellowstone -Snake River Plain basalts. From east to west, 
Yellowstone, Snake River Plain, Owyhee Plateau and the Oregon 
High Lava Plains (Fig. 3) are thought to reflect a continental hotspot 
track formed by the Yellowstone plume'*^. Snake River Plain lavas 
have high ^He/^He (Fig. 3), typical of OIBs, which were used to 
support a mantle plume origin for these lavas^°"^^. We note, however, 
that 'He/^He and 207pb*/206p^>i. increase and ^'^^Fhr'^Fh decreases 
from west to east along this hotspot track^^ (Fig. 3), indicating that 
^He/^He is positively correlated with ^^^^^^ p^^^p})^ and negatively 
correlated with "°'Pb/'°'Pb, just as in the global OIB-CFB dataset 
(Fig. 2). Snake River Plain lavas with high 'He/'He, ^^'Ph'^nh'^ 
and low 206p|^/204p|^ were emplaced through Precambrian North 
American Craton, while the High Lava Plains lavas with low 
'He/^He, ^^^Ph'^ r^Ph'^ and high ^^^Pb/'^^Pb were emplaced through 
Mesozoic accreted terrains. Although these spatio-temporal trends 
have been interpreted to represent changes in the dynamics of 
plume-lithosphere interactions^""^^, the close adherence of these iso- 
topic patterns to basement geology and age suggests that the Pb and 
He isotopic systematics may in fact be derived from the North 
American lithosphere and that the isotopic variations reflect varia- 
tions in the age of the Pb and He sources. A simple way of generating 
such systematics is to tap an ancient sulfide-bearing reservoir in the 
cratonal part of the Yellowstone-Snake River Plain volcanic series. 

In summary, the negative ^He/^He-^°^Pb/^°^Pb and positive 
'He/^He-"°'Pb*/'''Pb* trends in global OIB-CFBs and in continental 
intraplate magmas require a reservoir characterized by low time- 
integrated U/Pb and U/He, resulting in unradiogenic Pb and He 
isotopic signatures. Unradiogenic He is traditionally interpreted to 
represent an undegassed and hence primitive, unprocessed reservoir. 
Here, we suggest that the unradiogenic Pb (low 206p|^/204p|^ j^^gj^ 
^°^Pb*/^°^Pb*) and He isotopic signatures represent ancient mantle 
signatures frozen in the form of sulfide liquids co-precipitated with 
large masses of mafic cumulates (pyroxenites) during major melting 
events associated with the formation of continents. In this regard, the 
^He/^He isotopic signature of the mantle may reflect the integrated 
effect of major crust-forming events, as argued by Parman^^ but for 
different reasons. The reason why most of the upper mantle and 
continental crust is more radiogenic than the bulk Earth is thus 
because there exists a large reservoir of maflc, sulfide-bearing cumu- 
lates now stored in the mantle or deep lithosphere and complement- 
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ary to the continental crust today. Most recently, it has been sug- 
gested that amphiboles might serve as potential reservoirs of He and 
other noble gases^^. While this may be true, the stability field of 
amphiboles in the mantle is rather limited and, without sulfide, there 



is no obvious connection with Pb. We conclude that sulfide-bearing 
lithologies provide the most internally consistent explanation for the 
coupled Pb-He isotopic systematics of major mantle and crustal 
reservoirs. In particular, Pb and He isotopic systematics of the mantle 
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are downloaded from GeoRoc database, and are in Table S3. 

reflect the integrated production and destruction of continents rather 
than the preservation of primordial, unprocessed parcels of mantle 
that have survived unscathed throughout Earth's history. It is thus 
worth exploring the role of sulfides during silicate differentiation of 
the Earth, and in this regard, it maybe worth examining the impact of 
sulfides on the geochemical signatures of mantle- derived magmas 
through time. However, because S contents in magmas are controlled 
by sulfide solubility and the formation of cumulates occurs at high 
temperatures, the impact of sulfides is not likely to be easily studied 
by S content or S isotopes in magmas. A more fruitful approach 
would be to investigate the secular evolution of chalcophile elements, 
such as Cu, Ag, Pt and Pd, in the mantle as the behaviors of these 
elements during melting are strongly influenced by the presence of 
sulfides. 
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